Effect of disorder on a graphene p-n junction 
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We propose the theory of transport in a gate-tunable graphene p-n junction, in which the gradient 
of the carrier density is controlled by the gate voltage. Depending on this gradient and on the 
density of charged impurities, the junction resistance is dominated by either diffusive or ballistic 
contribution. We find the conditions for observing ballistic transport and show that in existing 
devices they are satisfied only marginally. We also simulate numerically the trajectories of charge 
carriers and illustrate challenges in realizing more delicate ballistic effects, such as Veselago lensing. 

PACS numbers: 81.05.Uw, 73.63.-b, 73.40.Lq 

I. INTRODUCTION AND MAIN RESULTS 
A. Definition of the model 

Graphene is a new material whose unique eleetronic 
structure endows it with many unusual properties.^ A 
monolayer graphene is a gapless two-dimensional (2D) 
semiconductor with a massless electron-hole symmet- 
ric spectrum near the corners of the Brillouin zone, 
e(k) — ±.hv\k.\, where v « lO^cm/s. The concentration 
of these "Dirac" quasiparticles can be accurately con- 
trolled by the electric field effect. An exciting exper- 
imental development is the ability to apply such fields 
locally, by means of submicron gates. Using this tech- 
nique, graphene p-n junctions (GPNJ) have been recently 
demonstrated . "^^^ 

Within idealized treatments that neglect disorder 
and electron interactions, GPNJ were predicted to dis- 
play a number of intriguing phenomena. They in- 
clude Klein tunneling, ^"^^ Veselago lensing, microwavc- 
induced^^'^3 and Andreev^'''^^ reflection, as well as strong 
ballistic magnetoresistance.^'^^ Both quantitative and 
qualitative changes to these phenomena are expected 
when interactions and disorder are included in the model. 
For example, long-range Coulomb interactions lead to 
nonlinear screening in GPNJ, which can modify its re- 
sistance substantially.^^ The purpose of this paper is to 
investigate how the junction resistance is affected by dis- 
order. We show that in existing GPNJ this effect is in- 
deed strong and suggest what can be done to reduce it. 

We consider a generic model of an electrostatic GPNJ, 
in which a grounded graphene sheet in the x-y plane 
is controlled by two coplanar metallic gates with volt- 
ages Vi and V2. The gates are separated by distance b 
from graphene and a distance 2d from each other. Un- 
der a symmetric gate bias, V2 — ~Vi — V (Fig. 1), the 
graphene carrier density n{x) varies linearly in the mid- 
dle of the junction (x = 0), 

n{x) ~ n' X , \x\ D = max {6, d} , (1) 

and tends to its limiting values ±no a-t |a;| 3> D. Here n' 



FIG. f : (a) Device geometry. Graphene (thin line) lies in the 
2 = plane. The gates (thick lines) are in the z = b plane, 
(b) Electron density profile for d = Q.77b and the symmetric 
gate bias V2 ~ —Vi — V. (c) Same for d = 6.00 b. 



is the density gradient^^ at a; = 0. 

Our assumptions about disorder require a brief discus- 
sion. At present, the nature of disorder in graphene is 
not completely understood.^ Our knowledge of it derives 
mainly from the measurements of the transport mobil- 
ity /i. For a sample with a macroscopically homogeneous 
carrier concentration n and resistivity p, the mobility is 
defined by 



1 



e\n\p{7i) 



(2) 



A remarkable fact that holds true for nearly all experi- 
ments on graphene is that /i(n) is observed to be approx- 
imately constant away from the charge-neutrality point 
n — 0. Rather than entering a debate on the microscopic 
origin of this behavior, we adopt it on phenomenologi- 
cal grounds. We can do so because the derivation below 
applies regardless of the exact microscopic origin of the 
constant mobility. 

It is convenient to define parameter rii of dimension of 
concentration by 



e 

hiJ, 



const , 



(3) 
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then the resistivity p(n) can be written as 

/ \ ^ I I 



(4) 



Below we will also need the carrier mean free path I, 
which is related to the conductivity in a standard way: 
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{2kFl): 



(5) 



where kpin) = \J 7r|n| 
Eq. (4), we find 

l{n) = 



is the Fermi wave vector. Using 



I'Krii 



\n\ > Hi . 



(6) 



The inequality \n\ » rii in Eqs. (4) and (6) is stipulated 

by another phenomcnological observation: the satura- 
tion of p{n) at a finite value p 

max 

h/e^ at low carrier 

densities. ^'^^ 

As wc mentioned, our main results can be obtained 
without knowing the microscopic origin of Eq. (4) and 
Pmax- Nevertheless, it is useful to have in mind a concrete 
model that may clarify the physical meaning of parame- 
ter m. One such actively discussed model assumes that 
the mobility is limited by charged impurities located in a 
close proximity to the graphene sheet. ^"'^^ An impurity 
of a unit charge acts as a scatterer with the transport 
cross-section^°~^^ 



A = 2'jTC2{a)/kF , 



(7) 



where C2 = 7ra^/2 for a ^ 1 (graphene on large-K sub- 
strate) and C2 ~ 0.1 for a 1 (Si02 substrate).^'' Here 
a = e^/nhv is the dimensionless strength of Coulomb in- 
teractions and K is the effective dielectric constant. If the 
charged impurities have an average surface concentration 
Ni, then I = l/(A^jA). Comparing this with Eq. (6), one 
indeed arrives at Eq. (3) with 



Hi = C2{a)Ni . 



(8) 



This argument has a considerable appeal and is sup- 
ported by recent experiments.^^ 



B. Results 

To isolate the transport properties specific to GPNJ 
we follow the procedure introduced by experimentalists,^ 
and compute the difference of the total resistance i?tot of 
the device in the p-n [Fig. 1(a)] and the n-n states: 



R = i?tot|y2=- 



Vi = V 



R 



.tot|y2=Vi=i^ • 



(9) 



This allows to largely eliminate the contribution of the 
bulk regions |a:| > D. Our results are then as follows. 
We find two qualitatively different regimes, depending 
on magnitude of the dimensionless parameter 



3/2 



(10) 



For small /3 (high disorder or low density gradient), the 
transport is purely diffusive, and the resistance of the 
GPNJ is given by 



/3<1: i?~2^ 



\n'\W 



■2/3 
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(11) 



where W is the width of the device in the y-direction and 
7 is defined by 



7 = |nf/3£)>l. 



(12) 



The condition 7 ^ 1, which is usually satisfied in 
experiment,**"^ ensures that the density n{x) varies 
across the GPNJ slowly enough, D = max{rf, b} » 
^^"'^('^^0)7 to justify its evaluation by means of classical 
electrostatics.-^"^ Equation (11) is written for /3^/^7 ^ 1, 
i.e., for no » rij, when the GPNJ is still well-defined de- 
spite random fluctuations of the electron density n{x, y) 
due to disorder. 

In the opposite regime (large (3 or low disorder) the 
GPNJ resistance 

/3 > 1 : R = Rbed + Rdii (13a) 
is the sum of the ballistic and the diffusive contributions, 
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(13b) 
(13c) 



Equations (11) and (13c) are valid with logarithmic ac- 
curacy^^ and match at /3 3. The ballistic contribution 
dominates, R ~ i?bai ^ -^dif > provided 



2ai/*5 , / 47r7 
m ' 



Cl 



3/2 



(14) 



Realistically, the logarithmically "large" threshold /3* 
here is about 10. In recent experiments"'''^ /? is of the 
same order of magnitude. So, they are presumably in 
the crossover region iibai ~ -Rdif- To move deeper into 
the ballistic regime one needs cither a larger concentra- 
tion gradient \n'\ or a higher mobility /i. 

The rest of the paper is divided into three sections. 
In Sec. II wc give the analytical derivation of the above 
results. In Sec. Ill we illustrate them by numerical simu- 
lations. Finally, in Sec. IV we discuss their implications 
for ongoing experimental work. 



II. DERIVATION 

This section is organized as follows. First, we con- 
sider electrostatics of the gate-tunable junction. Next, 
we study separately the ballistic and the diffusive contri- 
butions to the transport. Finally, we combine them to 
arrive at total expression for the resistance of a GPNJ. 
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density gradient sharply increases near the gate edges. In 
those regions, n{x) is dictated by the nearest gate (similar 
to the case studied in Ref. 17), and one can show that^* 



dn 
dx 



27eP 



max{|yi|,|y2|}, (21) 



FIG. 2: A sketch of the characteristic lengthscales in a GPNJ 
for the limiting cases of small and large f3. Only the a; > 
side of the junction is shown. The diffusive region is hatched. 
Parameters Xtun and are indicated by the dashed lines in 
the regimes /3 <g 1 and /3 ^ 1 respectively, as they do not 
have direct physical meaning in these cases. 



A. Electrostatics 

Electron density in graphene is related to the elec- 
trostatic potential ^{x, z) by the Gauss law, n{x) = 
{K/4:TTe)dz^{x, +0). To find <i> and n we can treat 
graphene as an ideal conductor. (For the discussion of 
this approximation, see Refs. 17 and 27.) The calcula- 
tion can be done using the conformal mapping 



2w + In ( = - {x + iz). 

a — w J b 



(15) 



which transforms the upper half-plane 2 > with the 
branch cuts along the gates [cf. Fig. 1(a)] to the upper 
half- plane of a complex variable w = w{x, z). Here a is 
found from 



y/a{a + 1) + ln(Va + Va + 1 ) = 7rd/(26) . 



(16) 



The graphene sheet, the left gate, and the right gate 
are mapped onto the intervals —a < w < a, w < —a, 
and w > a, respectively, of the real axis. Therefore, the 
sought potential is given by 

<^{x, z) = (I/tt) Im [V^i ln(a + w) - V2 ln(a - w)] . (17) 

Using these equations and simple algebra, we find 

^ {V2 + Vi)a + {V2-Vi)w{x) 



n(x) 



Sireb 



a{a + 1) — w^{x) 



(18) 



where w{x) stands for the real quantity w{x,z = 0) de- 
fined by Eq. (15). For the symmetric gate bias we obtain 



n{x) 



no{V)w{x) 



V2 



(19) 



a{a -I- 1) — w'^{x) 
in which case the density gradient at x = is given by 



TT no{V) 
26(l + a)2 



A'Keb 



(20) 



Two examples of n(x) computed according to Eqs. (15) 
and (19) are plotted in Fig. 1(b) and (c). In both cases 
the linear dependence n ~ n'x extends up to '-^ D. 
However, for widely separated gates [Fig. 1 (c)], the local 



B. Ballistic resistance 

The resistance i?bai of a clean GPNJ is related® to the 
electric field at the p-n interface. To compute this field 
one has to go beyond electrostatics of ideal conductors 
and take into account nonlinear screening at the p-n in- 
terface. Equation (13b) for i?bai was derived from this 
analysis in Ref. 17. In the case a ~ 1, the result for i?bai 
can be qualitatively understood as the ballistic resistance 
of a system with Wkp{n{xtna)) transmitting channels: 



1 



h Xtu 



Here the effective "width" of the p-n interface 

^tun 



a-i/6|„'|-i/3 



(22) 



(23) 



is found from the condition that it is of the order of the 
quantum uncertainty in the quasiparticlc coordinate, 



Xt^, 



n{xtun)) 



(24) 



(In Ref. 17, xtun was denoted by xtf ) The quasiparticles 
that manage to get inside the strip |x| < Xtun cross the 
p-n boundary without tunneling suppression.® 

Below we consider the resistance (9) of a symmetrically 
biased GPNJ, V2 = —Vi = V. The transport is either 
diffusive or ballistic depending on the gradient (10). 



C. Purely diffusive transport, /3 <C 1 

The derivation is based on treating p(n(a;)) as the lo- 
cal x-dependent resistivity. This is justified provided the 
concentration gradient is sufficiently small, such that 



l{n)\dxn\ < n . 



(25) 



Using Eq. (6), one can easily check that for (3^1 the 
condition (25) is satisfied at all |x| 3> Xfic where Xfic is 
defined by 



XRc 



nj\n'\ 



(26) 



see also Fig. 2(a). At such distances Eq. (4) is still 
valid. On the other hand, in the strip |a;| < xnc-, we 
have \n{x)\ < n^, so that Eq. (4) does not apply.^® Since 
the transport remains diffusive in the strip \x\ < xnc (cer- 
tainly, it cannot be ballistic because of strong disorder'^"), 
we can assume that the corresponding local resistivity is 
of the order of its bulk value Pmax ~ h/e^ at the charge 



4 



neutrality point. This allows us to estimate the resistance 
of this region as 



(27) 



According to our definition (9), the GPNJ resistance is 

the difference of the total resistances in the p-n and n-n 
configurations. It is convenient to write it as = T^-dif (0), 
where 

oo 

^dif(a;) = ^ jdx [p{x)\y^^_y^ - p(i)|y^=+yj . (28) 



Using Eqs. (15), (19), (4), and the expression 



n{x) 



no{V)a 



i{a + 1) — w'^{x) ' 



Vi=V2 = V, 



(29) 



for the charge profile in the n-n state that follows from 
Eq. (18), the integral in Eq. (28) can be transformed to 



Tldifix) ~ 2 — 



\W 



In 



no 



(a + l)|n'|a;_ 



> 



XRc 



(30) 

The total resistance is 7?.dif (a;flc) + Rrc, which leads to 
Eq. (11). Note that the effect of Ruc is only to modify 
the numerical factor in the argument of the logarithm in 
the final expression. For sufficiently long junctions, this 
logarithm is large, cf. Eq. (11), and so our crude esti- 
mate of Rfic is quite acceptable. This can be understood 
by realizing that in a long junction the resistance of the 
\x\ < Xfic strip is much smaller than that of the rest of 
the system. In shorter devices, the contribution of this 
"fluctuating strip" can be significant, and so a more ac- 
curate evaluation of TZdu iu Eq. (28) may be necessary. 
For example, one may want to perform the integration in 
Eq. (28) numerically using the experimentally measured 
dependence p{n) instead of Eq. (4). 



D. Co-existence of ballistic and diffusive transport, 

/3> 1 

Here the carrier density n{x) varies with x more 
rapidly. As a result, the diffusive approximation breaks 
down inside the strip |a;| < Xbai whose width is given by 
the condition ^[n(a;bai)] ~ aibai, i-e.. 



Airnj 



(31) 



The carrier density at x = a;bai is still high, n(xbai) ^ rii, 
so that at \x\ > Xha\ Eq. (4) applies. Thus, the diffu- 
sive contribution to the resistance is _Rdif — ^dif (■'z^bai), 
leading to Eq. (13c). [The extra factor /3~^ under the 
logarithm in Eq. (13c) vs. (11) comes from aJbai ^ f^'^xnc- 
Note, however, that xrc has no direct physical meaning 
if /3» 1.] 



In contrast, within the strip |.x| < .Tbai the transport is 
ballistic: the local mean-free path l[n{x)] nominally ex- 
ceeds \x\, so that quasiparticles reach the p-n interface 
largely without experiencing impurity scattering. We 
now note that the tunneling strip (31) is located deep 
inside this ballistic region. 



4n 



Q,l/6 134/3 



< a;bai 



(32) 



see also Fig. 2(b). Therefore, the transmission problem 
is reduced to the clean case,^'' yielding Eq. (13b) for the 
ballistic resistance i?bai- Due to the large logarithmic 
factor in i?dif [Eq. (13c)], the ballistic contribution in 
Eq. (13a) starts to dominate the diffusive one only when 
(3 exceeds a logarithmically large threshold /?*, Eq. (14). 



III. NUMERICAL SIMULATIONS 

In this section we illustrate and support the above an- 
alytical results by numerical simulations. In particular, 
we show that the criterion (3 ^ (3^, [Eq. (14)] guaran- 
tees only that the total resistance of the junction R is 
given by the formula derived for a disorder- free GPNJ, 
Eq. (13b). Realization of other ballistic phenomena may 
demand cleaner systems, see below. 

To get intuition above the nature of transport at 
/3 > /3* we studied semiclassical trajectories of the quasi- 
particles in a GPNJ by numerically solving the following 
relativistic equations of motion 



r = t;p/|p|, p = V|$(r)|. 



For illustrative purposes, we adopted the potential 

Q. 

3 



$(r) = — sgn {x)\/Tm'\x\ + 



(33) 



• (34) 



Here the first term models the potential induced by the 
gates"'^^ and the second term represents the potential cre- 
ated by impurity charges Qj = ±1 with coordinates 
{Tj,Zj). This expression assumes a = e^/K = Ti = v = l 
and neglects, for simplicity, the screening of these impu- 
rities by the electrons in graphene. We estimate that this 
entails C2 ^ 1 in Eq. (8). 

Other parameters of the simulation were as follows. 
The 2;-coordinates of all the impurities were set to Zj = 
0.01 in some arbitrary length units. The in-plane co- 
ordinates of the impurities were chosen randomly in- 
side the square jyl < 100 straddling the p-n in- 
terface. The total impurity number was 300, so that 
ni ~ 300/200^ = 0.0075. In the field of view < 60, 
\y\ < 80 of Fig. 3, 126 of these impurities are seen. The 
density gradient was set to be n' = 0.25, which makes 

parameter j3 quite large: j3 = 0.25/n^^^ ~ 400. 

In Fig. 3 we show 51 electron trajectories computed 

by standard numerical algorithms. ^"'^ The energy for all 
trajectories was fixed at zero and the starting point was 
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FIG. 3: Semiclassical trajectories of quasiparticles in a dis- 
ordered GPNJ with 3 S> 1. The trajectories start at the 
midpoint of the bottom edge, which belongs to the n-region. 
The p-region (shaded) occupies the upper half of the figure. 
The open (filled) circles arc in-plane Coulomb impurities of 
negative (positive) charge. The trajectories that carry cur- 
rent across the p-n interface are shown by dark lines on the 
n-sido and tapering white lines on the p-side. The variable 
width of these lines is the local Fermi wavelength 2nh/p{v). 
The thin lines are examples of the trajectories reflected from 
the interface. 



set to X = xo = —60, y = 0. The polar angles of the 
initial velocities formed an equidistant set spanning the 
interval (—■ 7r/5, 7r/5). 

From Eq. (31) we estimate .Xbai ^ 350. Therefore, the 
injection point is deep inside the ballistic strip, Xq <C 
Xbai- Simultaneously, xq ^ Xtun ~ 2, cf. Eq. (23), so 
that the semiclassical approximation (33) is legitimate. 

As evident from Fig. 3, the average distance between 
collisions of quasiparticles with impurities exceeds the 
distance xq from the injection point to the interface. 
Thus, in agreement with the above estimates, electrons 
can propagate across the interface according to the for- 
mulas derived for the disorder- free system. A closely 
related observation is that for the chosen parameters 
there are many points along the interface not "blocked" 
by the impurities. 

On the other hand, even for such large /? there is no 
evidence for the recently proposed^^ Vesclago lensing ef- 
fect: a self- focusing of holes into a point (xo, 0), a mirror 
image of the injection spot (— xo,0). The primary diffi- 
culties with observing this focusing effect are apparently 
as follows. First, even in the absence of any disorder, only 
a small fraction of electrons can penetrate through the 
p-n interface: the transverse momenta of such electrons 
must satisfy the condition^'^'' 



< 



^tun 



(35) 



Such momenta are much smaller than the typical ones 

h 



Py 



n'xo 



(36) 



Furthermore, scattering of an electron by a Coulomb 
impurity typically deflects the electron trajectory by a 
substantial angle. Therefore, the lensing additionally 
requires that a narrow fan of trajectories defined by 
Eq. (35) does not undergo impurity scattering. The 
width of this fan in real space is ^ ^/xoXtun ■ Therefore, 
the condition on xq becomes 



niXo^/XoXtun ^ 1 • 



(37) 



Accordingly, the injection and collection contacts must 
be placed no further than the distance 



1 



47r 



2/3^1/3 



a^bal 



(38) 



from the interface, which may be considerably smaller 
than Xbai- Indeed, the absence of a discernible Veselago 
lensing in Fig. 3 is in agreement with our estimates: since 
xo = 60 and Xiens ~ 20 [cf. Eq. (38)], we are not yet in 
the regime xq <C a;iens- 



IV. DISCUSSION AND CONCLUSIONS 

In this final section of the paper we discuss geometrical 
requirements imposed by the criterion (14) in actual ex- 
periments. Using a reahstic number /x ~ 2, 500 cm?/ (V-s) 
in Eq. (3), we get ~ 1.0 x lO^^cm"^. Such Ui can 
be achieved if the transport mobility is limited by, e.g., 
charged impurities of concentration Ni ~ 10^^ cm~^ [as- 
suming C2 ~ 0.1 in Eq. (8)]. 

We consider first the case of a narrow gap between 
the gates, 6 « d [Fig. 1(b)], where a « 0.4. Taking 
a 1, 6 Ri 50 nm, and no ~ 2 x lO^^cm"^, similar to 
those of Ref. 5, for the above rii we get /? ^ 10. Some 
evidence for the ballistic transport was indeed seen under 
such conditions.^ On the other hand, observing Veselago 
lensing^^ seems rather challenging: it requires placing the 
injection and collection contacts within ~ 10 nm from 
each other, cf. Eq. (38). 

Next, in the case of widely separated gates, d = 1 /xm 
(and the same b = 50 nm) we get /3 w 0.1 even for a 



very high maximum density no = 10 cm^ . In order to 
observe ballistic transport in this device, the suggested 
setup should be somewhat modified. For example, using 
a backgate one can introduce a uniform offset of the elec- 
tron density n(x), which would shift the location of the 
p-n interface away from the x = point and closer to the 
edge of either one of the gates, as discussed in Ref. 17. In 
this manner the density gradient n' at the GPNJ can be 
ramped up to its maximum value (21), yielding f3 similar 
to that in a narrow-gap device. 

Although we considered a particular junction geome- 
try. Fig. 1(a), our treatment can be readily extended to 
characterize transmission in any GPNJ with a smoothly 
varying electron density, 7 » 1. The basic steps are as 
follows: (i) find the carrier density gradient n' at the p-n 
interface, (ii) compute /3 from Eq. (10), (iii) determine. 
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based on the criterion (14), whether the device is dif- 
fusive or baUistic, and finaUy (iv) find the diffusive and 
ballistic contributions from Eqs. (ll)-(13b). [Formula for 
i?dif can be further refined if the integration in Eq. (28) is 
done numerically using an accurately measured density 
dependence of the bulk resistivity p{n).] 

To conclude, disorder can strongly inhibit the ballis- 
tic transport regime in grapliene field-effect devices. In 
recent expcrimcnts^^^ on graphcne p-n junctions, this 
regime was reached only marginally at best. For bal- 
listic devices one should aim at larger electron density 
gradients n' and higher mobilities to satisfy the condi- 
tion (14). Note that if the primary source of disorder are 
charged impurities, then the requirement on n' becomes 
less stringent for substrates of high dielectric constant 



K. In this case, on the one hand, n' is larger for the 
same gate voltage and, on the other hand, the influence 
of Coulomb scattering is smaller, C2 oc oc 
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